Asynchronous and corrected-asynchronous finite difference solutions of PDEs on MIMD multiprocessors

A major problem in achieving significant speed-up on parallel machines is the overhead involved with synchronizing the concurrent processes. Removing the synchronization constraint has the potential of speeding up the computation, while maintaining greater computation flexibility (e.g. differences in processors speed; differences in the data input to processors). We construct asynchronous (AS) finite difference schemes for the solution of PDEs by removing the synchronization constraint. We analyze the numerical properties of these schemes. Based on the analysis, we develop corrected-asynchronous (CA) finite difference schemes which are specifically constructed for an asynchronous processing. We present asynchronous (AS) and corrected-asynchronous (CA) finite difference schemes for the multi-dimensional heat equation. Although our discussion concentrates on the Euler scheme it should serve only as a sample, as it can be extended to other schemes and other PDEs.These schemes are implemented on the shared-memory multi-userSequent Balance machine. Numerical results for one and two dimensional problems are presented. It is shown experimentally that synchronization penalty can be about 50% of run time: in most cases, the asynchronous scheme runs twice as fast as the parallel synchronous scheme. In general, the efficiency of the parallel schemes increases with processor load, with the time-level, and with the problem dimension. The efficiency of the AS may reach 90% and over, but it provides accurate results only for steady-state values. The CA, on the other hand, is less efficient but provides more accurate results for intermediate (non steady-state) values. The results show the potential of developing asynchronous finite deference schemes for steady-state as well as non steadystate problems.This research was partially supported by a grant from The Basic Research Foundation administrated by The Israel Academy of Sciences and Humanities.A reduced version of the paper was presented at the 4th SIAM Conference on Parallel Processing for Scientific Computing, Dec. 11–13, 1989, Chicago, USA.The work by this author was supported by research grant 337 of the Israeli National Council for Research and Development in the years 1990–1991.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NASI-18107 while the author was in residence at the Institute for Computer Applications in Sciences and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665, USA.     

Monte carlo studies of phase transitions in polymer blends and block copolymer melts

Summary The unmixing transition of both symmetrical polymer blends AB (i.e. chain lengthsN _A=N _B=N) and asymmetrical ones (N _B/N _A=2,3) is studied by large-scale Monte Carlo simulations of the bond fluctuation model. Combination of semi-grand-canonical simulation techniques, ?histogram reweighting? and finitesize scaling allows an accurate location of the coexistence curve in the critical region. The variation of the critical temperature with chain length (N) is studied and compared to theoretical predictions. For the symmetrical case, use of chain lengths up toN=512 allows a rough estimation of crossover scaling functions for the crossover from Ising to mean-field exponents. The order-disorder transitions in melts of both symmetric (compositionf=N _A/(N _A+N _B)=1/2) and asymmetric (f=3/4) block copolymers is studied for very short chains (16≤N≤60). The interplay between structure and chain configuration is emphasized. Qualitative evidence for ?dumbell formation? of chains and vacancy enrichment in A-B-interfaces and near hard walls is presented. Paper presented at the I International Conference on Scaling Concepts and Complex Fluids, Copanello, Italy, July 4–8, 1994.     

非齐次Toda晶格的对称,精确解和可积性

研究非齐次Toda晶格，即一类非齐次非线性微分差分方程的对称与可积性。给出了这一类方程的Lie点对称，条件对称和精确解。给出这类方程与Toda晶格之间的可逆点变换，从而表明这一类方程是可积的。     

N维最弱受约束电子势模型的精确解

精确地求解了N维最弱受约束电子势模型束缚态的Schrodinger方程,得到了其能级和相对应的归一化径向波函数.在此基础上,运用广义拉盖尔多项式的级数展开得到了N维最弱受约束电子势模型径向矩阵元的通项表达式.     

Caloric properties of dilute nonelectrolyte solutions: A success story

Selected aspects of the thermodynamics of very dilute solutions of gases in liquids, in particular aqueous solutions, are reviewed and connected with recent high-precision experimental techniques [vapor-liquid equilibrium measurements (VLE), calorimetry and densimetry]. Some of the problems encountered in data reduction and data correlation over large temperature ranges, including the critical region, are discussed. The focus is on caloric properties, such as partial molar enthalpy changes on solution, ΔH^∞₂, and partial molar heat capacity changes on solution, ΔC^∞_P,2: direct calorimetric methods are compared with indirect methods based on VLE studies as a function of temperature (van't Hoff approach).     

具奇异边界条件的一类拟线性椭圆型方程

得到了一类拟线性椭圆型方程爆破解的存在性。     

NON-NEGATIVE SOLUTIONS FOR A CLASS OF NON-POSITONE PROBLEMS

1IntroduInthispaperwestudyill(.)(l}xlsl(.)ll(:()noTI-1)ega,t]v(3s(>1lltionsfort,h(i1,woilo]lit,boundaryvalueproblemsu\ry;l;(:i.t>(v,>\,:t),:j:=fti-i'...ifn(j)>l),A>0isdcollstanta,f,dj'e(!'l,.,ti,fi<,,f(())<0(l.:i)111.of)i(lfllsifs{.1.l:)-(.I.2)havel,ecnsi…     

The new identity for the scattering matrix of exactly solvable models

We discovered a simple quadratic equation, which relates scattering phases of particles on Fermi surface. We consider one-dimensional Bose gas and XXZ Heisenberg quantum spin chain. Received: 4 December 1997 / Accepted: 17 March 1998     

On the Quasiconcave Bilevel Programming Problem

Bilevel programming involves two optimization problems where the constraint region of the first-level problem is implicitly determined by another optimization problem. In this paper, we consider the case in which both objective functions are quasiconcave and the constraint region common to both levels is a polyhedron. First, it is proved that this problem is equivalent to minimizing a quasiconcave function over a feasible region comprised of connected faces of the polyhedron. Consequently, there is an extreme point of the polyhedron that solves the problem. Finally, it is shown that this model includes the most important case where the objective functions are ratios of concave and convex functions     

Mathematical modelling of rock bolt systems I

The main goal of the paper is to give a variational formulation of the behaviour of bolt systems in rock mass. The problem arises in geomechanics where bolt systems are applied to reinforce underground openings by inserting steel bars or cables. After giving a variational formulation, we prove the existence and uniqueness and some other properties.